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La presente invention concerne une methode de modelisation permettant de 
predire en fonction du temps la composition detail lee de fluides produits par un 
gisement souterrain en cours de production, et notamment d'un gisement 
5 d'hydrocarbures en cours de production. 

Un tel modele ~permet a i'ingenieur de gisement d'&ablir des profils 
compositionnels detailles qui, sont necessaires en particulier au dimensionnement 
pr6visionnel et a la gestion d'installations de surface, telles que separateurs, unites de 
traitement, lignes de transport, etc., done utiles en ing^nierie de surface et de 
10 procedes. 

Etat de la technique 

La modelisation des ecoulements dans un reservoir pettolier est fondee 
essentiellement sur Implication au reservoir prealablement maille* (ou a une portion 
de celui-ci) de la loi de Darcy bien connue decrivant Tecoulement de fluides en 
15 milieux poreux, de la loi de conservation des masses dans chaque unite* de volume, 
des relations thermodynamiques regissant Involution des parametres de phase des 
fluides tels que la viscosite, la masse volumique, sur des conditions initiales, sur des 
conditions aux li mites de fermeture de la structure, et sur des conditions aux puits. 

Le modele dit « Black Oil » designe ci-apres en abrege* par B.O., est Tun des 
20 plus utilises en simulation petroliere. II permet de decrire un ecoulement 
tridimensionnel et triphasique (eau-huile-gaz) compressible. Les effluents petroliers 
intervenant dans ce modele sont d^crits generalement par un constituant eau, et deux 
constituants pour le fluide de gisement, le mot de constituant recouvrant ici la notion 
de composant (comme H20 pour 1'eau) et la notion de pseudo-composant 
25 (regroupement de composants). Les constituants intervenant dans ce modele sont au 
nombre de trois : un constituant eau (E), un constituant hydrocarbure lourd (L), et un 
constituant hydrocarbure leger (V). Dans un modele de type B.O. dit « strict », le 



constituant (E) n'est present que dans la phase eau, le constituant (L) n'est present 
que dans la phase liquide hydrocarbure (appelee huile ou condensat), le constituant 
(V) se partage entre les phases hydrocarbures liquide et vapeur (phase appelee gaz). 
Meme si l'usage des modeles B.O. est a d^conseiller dans certains cas de gaz a 
5 condensats, il n'en reste pas moins applicable a un grand nombre de cas industriels. 

Un autre modele de simulation bien connu dit « compositionnel » est aussi 
utilise, ou les fluides hydrocarbures sont represented par un nombre plus elev6 de 
composants, au moins trois, souvent davantage. La modelisation des flux de ces 
fluides mieux d6taill6s se traduit par des temps de calcul tres importants (beaucoup 
10 plus longs que ceux requis pour une moderation de type B.O.) du fait du plus grand 
nombre de constituants mais aussi parce qu'il est souvent necessaire de require la 
taille des mailles pour limiter l'erreur numerique et par consequent d'augmenter le 
nombre de mailles. 

Pour des raisons pratiques, les fluides en place sont decrits comme constitues 
15 d'un nombre de composants ou pseudo-composants beaucoup plus r6duit que le 
nombre r£el de composants, de facon que les calculs de mod^lisauon puissent etre 
menes dans un delai raisonnable. Une composition r6duite a quelque 5 ou 10 pseudo- 
composants s'avere en general suffisante pour representer le comportement des 
fluides dans le gisement. 

20 Par la demande de brevet WO 99/42937 et aussi le papier de C. Leibovici et 

J. Barker "A Method for Delumping the Results of a Compositional Reservoir 
Simulation" SPE 49068, present^ a la SPE Annual Technical Conference and 
Exhibition New Orleans, 27-30 Septembre 1998, on connait une m&hode de 
prevision de 1'evolution de la composition detail lee au cours du temps a partir des 

25 calculs effectues dans une simulation de type compositionnel de fluides decrits par 
un certain nombre reduit de pseudo-composants (principe de representation 'groupee* 
ou 'regroupeV obtenue par une operation dite de 'lumping'), le nombre de 
composants etant au moins de trois. La methode permet done de prddire quels 


auraient ete les resultats d'une simulation de reservoir si Ton avait utilise" une 
modelisation finement detaillee ou les fluides sont represented par un nombre plus 
£leve de composants. C'est une operation connue des gens de 1'art sous V appellation 
de "delumping", que Ton designera par la suite aussi par "eclatement". 

Le principe de l'etape d'dclatement ant£rieure est de calculer le coefficient 
AD 0 et les n coefficients AD P (soit n+1 coefficients, n etant le nombre de parametres 
de l'equation d'6tat) d'une Equation generate connue, publi£e anterieurement dans un 
papier de C.F. Leibovici, E.H. Stenby, K. Knudsen, "A Consistent Procedure for 
Pseudo-Component Delumping", Fluid Phase Equilibria, 1996, 117, 225-232. : 

n 

(1) Ln(kj) = AD 0 + £aD p n pi 

P =i 

ou les Ilpi sont des parametres fixes de caracterisation du constituant i dans 
l'equation d'etat pour une representation thermodynamique donnee, a partir des 
constantes d'equilibre ki de chaque constituant de la representation thermodynamique 
regrouped calculees pendant la simulation compositionnelle dans chaque maille et a 
chaque pas de temps. Si est le nombre de composants de la representation 
thermodynamique regroupee, on dispose done de N rg equations pour determiner n+1 
coefficients. Une condition necessaire est done que N rg soit au moins 6gal a n+L 
Pour l'equation d'etat de Peng-Robinson a deux parametres, il faut done une 
representation thermodynamique regroupee a au moins trois composants. 

Une fois le coefficient AD 0 et les n coefficients AD P calcules, on les utilise 
pour calculer les constantes d'equilibre des composants de la representation 
thermodynamique detaillee (N rd composants) en appliquant l'equation (1) aux N„j 
composants avec leurs propres parametres fixes de caracterisation dans la 
representation thermodynamique detaillee. 

En utilisant a) les constantes d'equilibre pour la representation 
thermodynamique detaillee ainsi determinees, b) les flux entre chaque maille et aux 



puits, c) la fraction vapeur dans chaque maille issus de la simulation 
compositionnelle regrouped, et, d) la composition detaillee globale dans chaque 
maille et aux puits injecteurs en debut de pas de temps, on estime alors dans chaque 
maille la composition detaill6e de chaque phase hydrocarbure au pas de temps t et la 
5 composition detaillee globale de chaque maille au pas de temps (t+1) suivant. 

L'un des interets de cette methode est qu'il n'est pas n^cessaire, dans l'&ape 
d'eclatement (« delumping ») a chaque pas de temps, de resoudre l'equation d'6tat, 
que ce soit pour la representation regroupee ou la representation detaillee, ce qui 
constitue une 6conomie de temps de calcul. Un inconvenient de la methode est 
10 qu'elle n'est pas applicable aux simulations de type B.O, puisque les equations d'etat 
pratiques ont au moins deux parametres. 

Par la demande de brevet WO 98/5710726, on connait £galement une 
methode pour predire Involution de la composition detaillee au cours du temps a 
parti r des calcul s de flux effectues dans une simulation de type B.O. ou les phases 
15 hydrocarbures sont decrites par seulement deux composants (L) et (V). 
L'inconvenient de la methode est qu'elle requiert l'utilisation de l'equation d'etat de la 
representation detaillee a chaque pas de temps, et done d'etre couteuse en temps de 
calcul pendant l'etape de "delumping". 

La methode selon I 'invention 

20 La methode selon 1' invention combine une simulation previsionnelle du 

comportement du gisement par une modelisation de type B.O. avantageuse 
puisqu'elle ne n^cessite aucun recours a une equation d'6tat, et une 6tape 
d'eclatement ou de "delumping" qui ne necessite a aucun pas de temps une 
quelconque resolution d'une equation d'etat. 

25 Dans une modelisation de type B.O. on introduit, entre autres parametres, les 

parametres thermodynamiques des fluides (tels que viscosite, facteur de volume, 
densite, rapport de dissolution gaz-huile, etc.) dans un fichier des donndes d'entr6e 
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du simulateur (sous forme de tables, et ou par correlation, comme fonction de la 
pression, de la temperature si celle-ci varie) et, le cas ech£ant (pour le modele de 
type B.O. etendu), d'un parametre additionnel conservant une memoire de la 
composition du gaz tel par exemple que la densite du gaz). 

5 L'idee essentielle est ici d'introduire simplement dans les donnees d'entr£e les 

(n+1) colonnes ou correlations supplementaires correspondant aux differences AD P , 
par exemple en les incluant dans un fichier annexe au fichier principal utilise" par la 
simulation. II est 6galement possible d'introduire 2 x (n+1) colonnes ou correlations 
supplementaires correspondant aux parametres D p de chaque phase. Une condition 
10 pratique pour que des donnees puissent etre entrees sous forme de table est que les 
dependances fonctionnelles introduites soient suffisamment continues. 

Par la, on s'affranchit de la necessity du calcul des constantes d'equilibre par 
le modele. 

La methode de simulation selon V invention permet de predire en fonction du 
15 temps et dans au moins une zone thermodynamique, une composition detaillee d'un 
fluide produit par un reservoir et plus particulierement une composition detaillee 
d'un fluide contenu dans et produit par un gisement p&rolier dans lequel sont 
impl antes un ou plusieurs puits de production. 

Elle est caracterisee en ce que : 

20 - on discretise le gisement par un maillage, chaque maille contenant une ou 
plusieurs phases dont au moins une phase non aqueuse ; 

on determine la variation de parametres thermodynamiques des phases non 
aqueuses necessaires a une modelisation de type « Black Oil » (B.O.) au cours 
d'etapes d'un chemin thermodynamique suivi par les fluides dans le gisement 
25 (tels par exemple que la viscosite des phases, le rapport de dissolution du gaz 

dans l'huile, les facteurs de volume de l'huile ou de compressibilite du gaz, 
etc.) ; 


m 
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- on definit les fluides par une representation ddtaillee a N rd composants et/ou 
pseudo-composants ; 

- on determine au moins une equation d'etat a n parametres permettant de simuler 
au moins durant la phase de preparation des donnees d'entree (simulation qui est 

5 necessaire strictement que durant cette phase de preparation des donnees 

d'entree, et done non necessaire de . facon repetitive au cours ou apres la 
simulation de type B.O.), le comportement thermodynamique des fluides au 
cours des etapes d'un chemin thermodynamique suivi, permettant de generer, par 
zone ou domaine thermodynamique, n+1 relations fonctionnelles 
10 suppiementaires (sous forme de tables ou de correlations) ; 

- on traduit le comportement thermodynamique de chaque phase non aqueuse en 
des donnees d'entree convenant pour une representation thermodynamique de 
type B.O., donnees d'entree que Ton complete par les donnees d'entree 
suppiementaires convenant a l'operation d'eclatement ; 

15 - on realise une moderation de type B.O., permettant de determiner dans chaque 
maille et a des pas de temps consecutifs, des caracteristiques thermodynamiques 
de chaque phase non aqueuse et des donnees representatives de deplacement des 
phases dans le gisement ; et 

- on realise une operation d'eclatement pour obtenir la composition detailiee des 
20 fluides dans chacune des mailles. 

Pour caracteriser le chemin thermodynamique, on determine par exemple les 
variations de la composition des fluides et des proprietes des differentes phases. 

Au cours de la simulation de type B.O., on determine par exemple dans 
chaque maille des parametres thermodynamiques tels que la pression (pj ), la 
25 temperature (T7 ) si celle-ci varie, les saturations des phases hydrocarbure liquide 



(So^ ) et vapeur (Sg j ), les debits d'injection ou de production, et, pour chaque paire 
de mailles(j,h), les debits volumiques des phases liquide (Uo7h ) et vapeur (u^ h ). 

L' operation d'eclatement comporte par exemple une determination des 
constantes d'equilibre a partir de donnees d'entrees specifiques a 1'opexation 
5 d'eclatement, et une conversion de resultats exprim^s en volume par des n£sultats 
utilisables dans des equations de conservation des quantit^s mol aires ou massiques. 

Presentation sommaire des figures 

D'autres caracteristiques et avantages de la m£thode selon l'invention, 
apparaitront a la lecture de la description ci-apres d'un exemple non limitatif de 
10 realisation, en se reTerant aux dessins annexes ou : 

- les Fig. 1 a 3 montrent respectivement les variations en fonction de la pression de 
trois coefficients AD 0 , ADi, AD 2l au cours de la vaporisation diff£rentielle d'une 
huile de gisement simulee a l'aide de l'equation de Peng-Robinson a deux 
parametres. 

15 Description detaillee 

Les differentes etapes de la m^thode de modelisation permettant notamment 
de pr6dire en fonction du temps la composition detaillee de fluides produits par un 
gisement d'hydrocarbures en cours de production sont les suivantes : 

a) On represente tout d'abord le reservoir sous la forme d'un reseau de 
20 mailles (j) dont chacune constitue un volume elementaire rempli de fluide(s) sous 
forme de une ou plusieurs phases, avec au moins une phase non aqueuse. Les phases 
non aqueuses sont dites phases hydrocarbure meme si elles peuvent contenir certains 
composants autres que des hydrocarbures tel que l'azote, le dioxyde de carbone, 
Tanhydride sulfureux. 


8 

b) On etablit, d'une facon connue en soi, le comportement 
thermodynamique de la (des) phase(s) non aqueuse(s) necessaire a une 
representation thermodynamique de type B.O.. II s'agit de decrire la dependance de 
parametres de phase avec la pression, la temperature si celle varie en cours de 
5 simulation et, eventuellement, pour les moderations de type B.O. non strict, un 
indicateur devolution de la composition (par exemple le contenu en condensat 
dissous dans la phase gaz) au cours d'un chemin (Evolution) thermodynamique 
representatif de celui que va suivre le fluide dans le gisement. Le chemin 
thermodynamique le plus courant est celui de Toperation dite diffdrentielle au cours 

10 de laquelle la pression du melange est progressivement abaissee a temperature 
constante : sous le point de saturation, Tune des phases est evacuee a chaque palier 
de pression (pour la vaporisation differentielle d'une phase liquide, la phase gaz est 
evacuee); durant cette operation, les compositions et proprietes des phases 
(viscosit6s, masses volumiques, facteur de volume, rapport de dissolution gaz-huile 

15 de la phase liquide, etc.) sont determinees. 

II est a noter que Ton peut egalement traiter les cas qui requierent l'usage de 
plusieurs representations thermodynamiques « Black Oil » par exemple, si la 
composition en debut de simulation varie avec la profondeur ou lateralement a 
Tinterieur du reservoir, ou si Ton peut distinguer plusieurs chemins 
20 thermodynamiques locaux durant la moderation. Plusieurs zones ou domaines de 
variations des grandeurs thermodynamiques ou de composition, ddsignes sou vent par 
les specialistes comme zones thermodynamiques, pourront done etre definis et 
utilises. 

L'origine des donnees sur le comportement du fluide de gisement peut etre de 
25 nature expenmentale (experiences de laboratoire), elles peuvent etre issues de 
mesures sur champ, ou predites a partir de donnees frustres et d' hypotheses 
complementaires en simulant numenquement le comportement du fluide de 
gisement. Parmi les proprietes de phase ainsi habituellement decrites en entrees des 
modeiisations de type B.O. on trouve la viscosite des phases, le rapport de 
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dissolution du gaz dans l'huile, le facteur de volume de l'huile, le facteur de 
compressibilite du gaz. 

c) On deTinit, avant la simulation B.O. et pour chaque zone ou domaine 
thermodynamique, le fluide par une representation detaillee a composants et/ou 

5 pseudo-composants, 

d) Par zone thermodynamique, on d6termine une equation d'etat 
permettant de reproduire prealablement, avant la simulation B.O., le comportement 
thermodynamique du fluide au cours du chemin thermodynamique repr6sentatif de 
celui que va suivre le fluide dans le gisement. 

10 e) On reproduit avec l'equation d'etat les etapes (souvent des paliers de 

pression) du chemin thermodynamique de chaque zone, en particulier sous la 
pression de saturation ou deux phases sont presentes. Lors de cette simulation 
thermodynamique prealable, on utilise, a la pression de saturation du fluide et en 
dessous, a chaque etape du chemin thermodynamique, les constantes d*6quilibre 

15 obtenues a partir des r^sultats de la simulation, pour calculer le coefficient ADo et les 
n coefficients AD P de l'equation (1) correspondant a l'equation d'etat choisie pour 
decrire le comportement du fluide hydrocarbure par exemple en minimisant la 
fonction : 


(2) O(AD 0t AD 1 ,..,ADp,..,AD n )= X 


2 

Nrd 

AD 0 + £AD P npi - Ln(ki) 


i=l 


20 L'equation (1) ayant et^ obtenue a partir d'une equation sur la fugacite du 

constituant i dans une phase ou apparaissent les coefficients Do et les n coefficients 
Dp propres a la phase : 


r f. \ n 
(3) LnK =D 0 +EDpn pi) 


10 

on peut alternativement determiner ces coefficients qui apparaissent comme des 
parametres de phases, l'inconvenient etant d'augmenter le volume des donnees 
d'en trees. 

On dispose alors a chaque etape du chemin thermodynamique, privu comme 
representatif de celui que va suivre le fluide dans le gisement, de la d6pendance en 
fonction de la pression, de la temperature et, le cas echeant, d'autres indicateurs 
necessaires, des proprietes des phases necessaires a la simulation B.O., et, sous la 
pression de saturation, du coefficient ADo et des n coefficients AD P (ou des n+1 
parametres D 0 et D p , p variant de 1 a n, par phase) qui seront utilises pour calculer 
les constantes d'equilibre lors de l'etape de "delumping" decrite ci-apres. 

Les variations des differences AD P (p= 0 a 2) pour une operation 
differentielle sur une huile de gisement simulee a I'aide de Tequation de Peng- 
Robinson a deux parametres, sont illustrees aux Fig.l a 3. On peut constater que le 
comportement de ces parametres est suffisamment continu pour etre aisement 
15 introduit sous forme de table, ou alternativement par une correlation simple en 
fonction de la pression. Par exemple, dans le cas present pour ADo: 

ADo = 1.406231 ( 1 - LN(P/Psat) [ 1 + 0.038 x LN(P/Psat)]) 

Au stockage pour utilisation ulterieure de revolution des valeurs de ces 
coefficients en fonction des abscisses de pression, temperature, et indicateur 
20 compositionnel, il peut etre pratique d'ajouter 1'evolution des valeurs de masse 
volumique de Thuile dans les conditions thermodynamiques de chaque etape du 
chemin thermodynamique suivi, tant au-dessus de la pression de saturation que au- 
dessous. On peut choisir 6galement ou alternativement de stocker d'autres parametres 
des phases huile et/ou gaz, mais ce n'est pas indispensable. 

25 f) On traduit, d'une facon connue en soi, le comportement 

thermodynamique de la (des) phase(s) hydrocarbure(s) en une ou plusieurs 
representation(s) thermodynamique(s) de type B.O. dans les donnees d'entree de la 


5 


10 
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modelisation. Ces d^pendances sont decrites au-dessus de la pression de saturation 
du melange d'hydrocarbures et au-dessous de la pression de saturation, soit point par 
point dans des tables avec une methode d' interpolation et, le cas ech^ant, une 
methode d' extrapolation, soit par des correlations analytiques, soit par un melange 
5 de tables et de correlations. 

g) On realise, de facon connue en soi, une simulation de type B.O., ladite 
simulation permettant de calculer au moins dans chaque maille (j) et a des pas de 

m m 

temps consecutifs (m, m+1, etc.), la pression (pj ), la temperature (Tj ) (si celle-ci 

m m 

varie), les saturations des phases hydrocarbure liquide (Soj ) et vapeur (Sgj ), les 
10 debits d' injection ou de production, et, pour chaque paire de mailles (j,h), les debits 

m m 

volumiques des phases liquide (u Q jh ) et vapeur (ugj h ). 

h) On determine, a chaque pas de temps (m) et dans chaque maille (j), 
les constantes d'equilibre (kij ) du composant i a partir des dependances des 
coefficients (AD 0s ADi,..,AD p ,..,AD n ) de l'equation (1), ou des coefficients (D 0 , 

m 

15 Di,..,D p ,..,D n ) de l'equation (3) dcrite par phase, £tablies en e) avec la pression (pj ), 
la temperature (Tj 1 ) et avec, le cas £ch6ant en cas de modelisation B.O. non stricte 
un indicateur de composition comme le contenu en condensat dissous dans la phase 
gaz. 

m 

i) On determine, a chaque pas de temps (m), la fraction vaporisee (0j ) 
20 dans chaque maille (j), soit a partir des resultats de la simulation de type B.O. soit, 

pour une meilleure precision en cas de coexistence des phases huile et gaz, en 
resolvant l'equation de Rachford-Rice connue des specialistes qui est appliqu^e par 
exemple dans le document suivant : 

- Rachford, H.H.Jr et Rice, J.D.; "Procedure for Use of Electronic Digital 
25 Computers in Calculating Flash Vaporization Hydrocarbon Equilibrium", J. Pet. 

Technol., 1952, 14, 19, 
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a partir des fractions molaires de chaque composant i dans la composition detaillee 
globale (z™ ) du fluide hydrocarbure dans la maille (j) au pas de temps (m) 

Nrd 7 m (k m _ n 

j) On estime la composition detaillee de chaque phase hydrocarbure, a 
5 chaque pas de temps (m) et dans chaque maille (j), en utilisant les relations suivantes 


(5) xy = Z " m M « ~ lm pour la phase huile 

1 + (ley -1) 0j 

k™ z™ (k™ — 1) 

(6) vh = 11 2 11 m '' m pour la phase gaz 

v J y i 1 + (kij -l)9j 

k) On estime la density molaire de la phase huile £o7 dans la maille (j) 
10 au pas de temps (m) par exemple a partir de la masse volumique de l'huile poj , dans 
la maille (j) au P 35 ^ e temps (m), par la relation : 

m 

(7) £o7 = mmq" > MMo etant la masse molaire de la phase liquide 
hydrocarbure (elite 'huile'). 

Dans l'equation (7), la masse volumique po^ peut done etre issue des 
15 resultats au pas de temps (m) de la simulation B.O. ou obtenues a partir des donnees 
d'entree specifiques a l'etape d*£clatement stockees lors de la simulation 
thermodynamique pr^alable. MMo peut etre calcutee a partir des resultats de la 
simulation B.O. si la masse molaire du constituant hydrocarbure lourd est deTinie par 
defaut ou dans les donnees d'entrees utilisateur, ou a partir de : 
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(8) MMoj = X X U » ^ tant la masse molaire du composant ou 

i=l 

pseudo-composant i, et les quantity x™ ayant ete calculees dans l'6tape precedente. 

1) On estime la densite molaire de la phase gaz £g7 dans chaque maille 
(j) a chaque pas de temps (m), soit a parti r de donnees d'entree specifiques a Tetape 
5 d'eclatement stockees lors de la simulation thermodynamique pnSalable, soit a partir 

tn y m 

des resultats de la simulation B.O. en pression Pj , temperature Tj, par exemple 
avec, lorsque le facteur de compressibilite du gaz fait partie des donnees d'entnSe de 
la representation thermodynamique de type B.O. la relation suivante : 

m 

m P ■ m 

(9) = — hH — m , Zg i etant le facteur de compressibilite du gaz, R la 

Zgj R Tj J 

10 constante des gaz parfaits 

ou, par exemple, si c'est le facteur de volume du gaz qui fait partie des donnees 
d'entree de la representation thermodynamique de type B.O. avec la relation 
suivante : 

(10) £e? = — sr^ , Vmolst etant le volume molaire du gaz en conditions 

1 Bgj Vmolst 

15 standard et Bg™ etant le facteur de volume du gaz. 

m) On evalue dans chaque maille (j) a chaque pas de temps (m) les debits 
molaires des phases liquide (u c 7h ) et vapeur (u^L ) a partir des debits volumiques des 
phases liquide (i^h ) et vapeur (u^ h ) obtenus lors de la simulation B.O. et les densit^s 
molaires obtenues dans les etapes k) et 1). 

20 n) On estime au pas de temps m+1 la fraction molaire de chaque 

m+l 

composant ou pseudo composant i dans la composition detail lee globale (z y ) du 



fluide hydrocarbure dans la maille (j) connaissant la fraction molaire (zy ) au pas de 
temps m a parti r des equations ci-apres : 

mm m „m m m — m m m m 

zy N j - At (yy Sgj + xy Soj) - At 2, (yij' u sih + H' "ojn) 

he JO) 

m+1 

(11) ZU = tff 

m+1 m m in 4 ^ m m 

(12) N j = Nj - At (Sgj + Soj) - At £ (Ugjh + u ojh ) 

heJ(j) 

m 

5 l'equation (12) traduisant !e bilan molaire total hydrocarbure sur la maille Nj > 
etant la quantity molaire hydrocarbure contenue dans la maille (j) au pas de temps 
(m) en tenant compte des echanges de matiere avec toutes les mailles (h) voi sines de 
(j) qui forment l'ensemble J(j). 

Dans liquation (11), l'ecriture des termes (y™) et (x^) , - dans lesquels j'=j 
10 pour un debit de la maille (j) vers la maille (h) ou dans le puits, et j*=h pour un d6bit 
de la maille (h) vers la maille (j), et j' correspondant au fluide injecte dans le cas des 
puits d'injection, S etant alors ntSgatif suppose implicitement l'utilisation d'un 
schema simple amont pour les flux compositionnels. Une ecriture plus gendraie de 
ces termes est (y™ h ) et (x^) ou x^ et y™ h decrivent les compositions des phases 
15 liquide et gaz , obtenues d'une facon connue en soi, s'dcoulant entre les mailles (j) et 
(h). 

Les termes source (debits molaires) (S^) et (Soj) sont obtenus a partir des 
termes source (debits volumiques) (S^) et (S™) dans la simulation B.O. en 
reproduisant pour les termes source les Stapes j) a m). 

20 Les etapes h) a n) decrivent l'operation d'eclatement ou "delumping". Les 

details fournis sont bases sur des equations de conservation en moles ; on aurait pu 
tout aussi bien se baser sur des equations de conservation en masse. 
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Des lors que Ton sait decrire la composition detaillee dans chaque maille du 
modele de simulation de type B.O. a l'instant t, qui peut etre en particulier l'instant 
initial, on sait grace a la presente mdthode, decrire Involution de la composition 
detaillee dans chaque maille au cours du precede* d'exploitation mod^lisd dans un 
modele de type « Black Oil ». 
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REVENDICATIONS 

1) Methode de modelisation permettant de pr6dire en fonction du temps la 
composition detaillee des fluides produits par un gisement souterrain en cours de 
production, dans au moins une zone thermodynamique definie dans la modelisation, 
5 caracterisee en ce que : 

on discretise le gisement par un maillage, chaque maille contenant une ou 
plusieurs phases dont au moins une phase non aqueuse ; 

- on determine la variation de parametres thermodynamiques des phases non 
aqueuses necessaires a une modelisation de type « Black Oil » (B.O.) au cours 

10 d' Stapes d'un chemin thermodynamique suivi par les fluides dans le gisement ; 

- on deTinit les fluides par une representation detail lee a composants et/ou 
pseudo-composants ; 

on determine au moins une fonction d'etat a n parametres permettant de simuler 
au moins dans la phase de preparation des donnees d'entree, le comportement 
15 thermodynamique des fluides au cours des etapes du chemin thermodynamique 

suivi de facon a generer, dans chaque dite zone thermodynamique ou Ton 
souhaite effectuer une operation d'eclatement, des donnees d'entree 
suppiementaires ; 

- on traduit le comportement thermodynamique de chaque phase non aqueuse en 
20 des donnees d'entree convenant pour une representation thermodynamique de 

type B.O., donnees d'entree que Ton complete par les dites donnees d'entree 
suppiementaires convenant a une operation d'eclatement ; 

- on realise une modelisation de type B.O., permettant de determiner dans chaque 
maille et a des pas de temps consecutifs, des caract£ristiques thermodynamiques 

25 de chaque phase non aqueuse et des donnees representatives de deplacement des 

phases dans le gisement ; et 
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on realise une operation d'eclatement pour obtenir la composition detaillee des 
fluides dans chacune des mailles. 

2) Methode selon la re vendi cation 1, caracterisee en ce que Ton caracterise le 
chemin thermodynamique en determinant les variations de la composition des 
fluides et des proprietes des differentes phases. 

3) Methode selon la revendication 1 ou 2, caracterisee en ce que Ton 
determine dans chaque maille au cours de la simulation de type B.O., des parametres 
thermodynamiques tels que la pression (p™ ), la temperature (Tj 1 ), si celle-ci varie, 

m m 

les saturations des phases hydrocarbure liquide (Soj ) et vapeur (Sgj ), les debits 
d' injection ou de production, et, pour chaque paire de mailles(j,h), les debits 

volumiques des phases liquide (u 0 jh ) et vapeur (Ugjh). 

4) Methode selon Tune des revendications prec6dentes, caracterisee en ce que 
l'op6ration d'eclatement comporte une determination des constantes d'equilibre a 
parti r de donnees d'entr^es specifiques a l'operation d'eclatement, et une conversion 
de resultats exprimes en volume par des resultats utilisables dans des equations de 
conservation des quantites molaires ou massiques. 

5) Methode selon Tune des revendications precedentes, caracterisee en ce que 
Ton utilise chaque fonction d'etat pour generer n+1 relations fonctionnelles 
supplementaires sous forme de tables de donnees ou de correlations, que Ton inclut 
dans les donnees d' entree. 
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